Break— Junction Tunneling on MgB2 
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Tunneling data on magnesium diboride, MgB2, are reviewed with a particular focus on 
superconductor-insulator-superconductor (SIS) junctions formed by a break-junction method. The 
collective tunneling literature reveals two distinct energy scales, a large gap, Ai ^ 7.2 meV, close to 
the expected BCS value, and a small gap. As ~ 2.4 meV. The SIS break junctions show clearly that 
the small gap closes near the bulk critical temperature, Tc = 39 K. The SIS spectra allow proximity 
effects to be ruled out as the cause for the small gap and therefore make a strong case that MgB2 is 
a coupled, two-band superconductor. While the break junctions sometimes reveal parallel contribu- 
tions to the conductance from both bands, it is more often found that As dominates the spectra. In 
these cases, a subtle feature is observed near As + Al that is reminiscent of strong-coupling effects. 
This feature is consistent with quasiparticle scattering contributions to the interband coupling which 
provides an important insight into the nature of two-band superconductivity in MgB2. 

PACS numbers: 73.40.Gk, 74.50.+r, 74.70.Ad, 74.80.Fp 



INTRODUCTION 

Approximately one and a half years after the discovery 
101 of superconductivity in MgB2 a wealth of information 
about its properties has been collected. While its high 
critical temperature and simple crystal structure gener- 
ated an initial flurry of activity, there is now an increasing 
interest in the nature of superconductivity in MgB2, as 
it appears to be one of the rare examples of a two-band 
superconductor. 

Theoretically, such two-band behavior was proposed 
for MgB2 because the electronic system consists of two 
qualitatively different types of charge carriers, derived 
from boron tt- and cr-bands, respectively The tt- 
bands are three dimensional (3D), while the cr-bands are 
effectively restricted to two dimensions (2D). For conve- 
nience, these actual four bands are henceforward treated 
as two effective bands. Superconductivity is proposed 
to arise from electron-phonon coupling of the 2D band 
with a specific boron bond-stretching mode. That im- 
plies, that superconductivity originates from the cr-band 
and that superconductivity in the 7r-band is driven by 
that primary interaction. As a result, the properties of 
superconductivity are proposed to be different in both 
bands, provided the material is in the clean limit and the 
electronic systems do not strongly mix. The 2D band 
shows a large gap. A/;, ~ 7.2 meV, whereas the 3D band 
shows a small gap, A5 ~ 2.4 meV, both closing at a joint 
critical temperature, Tc- In the dirty limit, only one gap 
of intermediate magnitude is expected to be observed, 
closing at a reduced Tc. 

There now is sufficient experimental evidence to 
strongly support such a two-band scenario, and there is 
a growing consensus in the community to this end. How- 
ever, a large fraction of these experiments involved tun- 
neling spectroscopy which is susceptible to surface imper- 



fections. Specifically, surface proximity effects can mimic 
the effects of two-band superconductivity. In this paper, 
we will present an overview of tunneling results, and will 
use this together with our own experimental data to ar- 
gue in favor of the two-band nature of superconductivity 
in MgB2. 



TECHNIQUES OF TUNNELING 
SPECTROSCOPY 



Tunneling spectroscopy traditionally presents one of 
the most direct probes of the superconducting energy 
gap In superconductor-insulator-normal metal (SIN) 
tunnel junctions, the conductance at low temperature is 
equivalent to the density of states (DOS) near the Fermi 
energy, Ep. At finite temperature, the conductance may 
be calculated from: 

where N is the superconducting DOS and cr^r the conduc- 
tance of the junction in the normal state. To account for 
experimentally observed broadening in excess of thermal 
smearing, a smearing parameter V usually is introduced 
into the BCS DOS Hi: 



N{E) = SR 



V(|£;|-ir)2-A2 



(2) 



A widely used realization of such SIN junctions is the 
scanning tunneling microscope (STM), where the barrier 
is given by a short length of vacuum that separates the 
STM tip from the sample. This technique combines an 
ideal insulating barrier with high spatial resolution. On 
the other hand, quite frequently a tip is placed into di- 
rect contact to the sample, a setup generally known as 
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FIG. 1: Simulated temperature dependence of (a) an SIN junction with Z = (purely metallic), (b) Z = 0.5 (intermediate), 
(c) Z oo (pure tunneling) and (d) an SIS tunnel junction. Parameters used are: Ao — 2.4 meV and Tc = j-^fe = l^.S K. 
Curves are vertically shifted for clarity, temperature increment is 1 K. To avoid singularities, F = 1 fieV was included (using 
the formalism described in Ref. § for (a) and (b)). 



point contacts. Such systems usually lack fine spatial 
resolution, but they allow for a significantly lower con- 
tact resistance and can be used to realize a continuous 
transition from a tunnel junction (insulating barrier) to 
a metallic contact (no barrier). Metallic contacts show 
enhanced current at bias voltages lower than the super- 
conducting energy gap due to Andreev reflections Q , and 
hence they also allow one to determine the energy gap. 
In the intermediate case, both tunneling and Andreev 
reflections contribute to the transport over the barrier, 
a case that was analytically described by Blonder, Tin- 
kham and Klapwijk (BTK) @. 

Finally, symmetrical SIS junctions display very sharp 
features at twice the energy gap, that are a result of the 
convolution of two BCS DOSs: 

I=^J N{E)N{E + eV)[f{E)- f{E + eV)]dE. (3) 

Such SIS junctions may be used to study the energy gap 
and a major advantage is their ability to trace the gap 
feature to highest temperatures, as the coherence peaks 
prove to be fairly insensitive against thermal smearing. 

To demonstrate this. Fig. |^ compares the temperature 
evolution of simulated spectra assuming a 2.4 meV BCS 
superconductor with insignificant smearing (F = 1 /ieV) , 
for (a) a purely metallic, (b) an intermediate and (c) 
a pure tunnel junction, as well as (d) a symmetric SIS 
tunnel junction. Not only do the coherence peaks re- 
main sharp up to highest temperatures in the SIS case, 
their position furthermore directly scales with the clos- 
ing of the gap. In comparison, the peak position in the 
SIN tunneling case increases as T approaches Tc — even 
though the gap actually closes. A close relation between 
the SIS coherence peak position and A(r) is preserved 



even in the presence of significant smearing. We will 
make use of this peculiarity to precisely measure A5 (T) . 

EXPERIMENT 

A tremendous amount of tunneling work has been per- 
formed on MgB2, and even though such a narrow view 
certainly does not mirror the full depth of information 
obtained from these studies, we for now want to focus 
only on the inferred energy gap values. Table | gives a 
compilation of these values, and it immediately becomes 
evident that most studies agree with the presence of two 
gaps around As ~ 2.0 — 2.8 meV and Al = 7.0— 7.5meV. 
This observation of two gaps, as well as their values, are 
in nice agreement with theoretical predictions ||^, ^ , 
however, some more effort is needed to positively con- 
firm two-band superconductivity as the origin of these 
features. 

Tunneling spectroscopy by nature is a surface probe 
and it is important to carefully verify that its results 
represent the bulk properties of the sample. Degraded 
surfaces were repeatedly suggested as an origin for small 
gap values |39), and there still are some lingering con- 
cerns as to whether multiple gaps in tunneling are "real" . 
To answer to that, more information than the mere gap 
values is needed, and tunneling can readily supply such 
information. 

The polycrystalline samples used in the present work 
were prepared as described in Ref. ]4l| |. Tunneling spec- 
tra were taken using a home-built point-contact appa- 
ratus and a Au tip Q. This resulted in both SIN and 
SIS junctions, however, henceforward we will focus on 
the SIS break junctions. 
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Experiment 



Energy Gap / [meV] 



Reference 



STM 



2.0 



A = 5 - 7 





A 

As = 


2.3 


Al 


= 7.1 




As = 


3.5 


Al 


= 7.5 




Aq = 


2.2 


Al 


= 6.9 






5.0 


Az 


= 8.0 


MgB2-In, Ag Paint 


As = 


2.6 






Ma-Ro-Pt 


Ac; - 


1.7 


Ar 


— 7 


IVTg-Ro-Nh 




3 


Ar 


= 7 


MgR2-Nb, Ptir 


Aq = 


2.8 


Ar 


= 9.8 




Aq — 


2.8 


Ar 


— 7.0 




Ac — 


2.45 


At 


— 7 n 


MgB2-Au 


As = 


2-3 


Al 


= 6- 


MgR2-Au 


As = 


2.3 


Al 


= 6.2 


MgB2-Au 




A = 3 


- 4 




MgR2-Au, Pb 


As = 


2.7 


Ai 


= 7.1 


MgR2-Au, Pt, PtIr, In, Ag Paint 


As = 


2.8 


Ai 


= 7.1 




As = 


1.7 - 2 








As = 


2.5 


Al 


= 7.6 




Point Contact 



|22| 
|23| 
|24| 

125,0 
|27| 

|28,§,|§ 

l"l 

|32| 
|33| 
|34| 

|35,g,p7t 



Rreak Junction 



(SQUID) 



As = 2 - 2.25 
As = 2.02 



8.5-9.5 



N/A 



|37, fe| 

B B 



1 



|44| 

|45j^ 



Planar Junction 



Step-Edge 
Ramp 
Nanobridge 
Trench 



MgB2/Ag/MgB2 
MgB2/MgO/MgB2 
(SQUID) 
(SQUID) 



Metal Masked Ion Damage 



N/A 
N/A 
N/A 
N/A 
N/A 



|47 
149, 5C 



1 51 
1 54 



52, 



Sandwich 



MgB2/Ag 
MgB2/Pb 

MgB2/Al/Al2 03/Nb 
MgB2/AlN/NbN 
MgB2/MgO,/Au, Ag 



As = 1.75 
As = 2.2 
As = 2.95 
As = 2.5 



Al = 7.3 
Al = 8.2 



|56| 
|57| 
|5S| 
|59| 
|60| 



TABLE I: Energy gaps of MgB2 inferred from different tunneling experiments in the literature. Experiments are listed by 
technique as well as they are sorted by research group. For each group, only the most recent findings are given. For reference, 
all tunneling studies are listed, including a few that did not determine any energy gap value. 



TEMPERATURE DEPENDENCE OF As 

The first unexpected result from tunneling spec- 
troscopy on MgB2 was the observation of a small energy 
gap, A5 — 2.0 meV [||, much smaller than the weak 
coupling limit would allow for = 39 K. The simplest 
scenario to result in a reduced gap value is a surface layer 
of reduced superconductivity with a gap anywhere in be- 
tween the bulk value and zero. This scenario can posi- 
tively be ruled out by studying the temperature depen- 
dence of the small gap, since such a reduced gap should 
still scale with the local, reduced Tc, i.e. it should obey 
> 3.52 0. This is illustrated in Fig. |, where 



schematically the temperature dependence of an intrin- 
sic small gap is compared with the expectation for As in 
a coupled two-band model. 

For this reason, As(T) was measured in a number of 
tunneling studies, and it was concluded to close at or 
near the bulk Tc of 39 K. However, the unusual com- 
bination of a small gap and high temperatures results 
in significant smearing in SIN junctions that renders the 
determination of As{T) more difficult. In cases where 
two gaps are observed simultaneously the exact analy- 
sis of high-temperature data is further complicated. SIS 
junctions promise to provide a solution to this problem, 
as the sharpness and position of the coherence peaks are 
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FIG. 2: Temperature dependence for an uncoupled small gap 
(heavy line, which was used to compute Fig. hf) as compared 
to a small, coupled gap (thin line after Ref. p^], dashed line 
after Ref. 

insensitive to thermal smearing. 

The evolution of SIS conductance spectra with tem- 
perature (Fig. ^) also allows a clear identification of a 
junction as being SIS type, whereas the low-temperature 
data might be inconclusive. To illustrate this, Fig. ^ 
shows a junction that originally was interpreted in terms 
of SIN tunneling using A = 4.1 meV and no smearing, 
but no temperature dependent data was taken |o|. 
This fit is compared to an SIS model using A = 2.3 meV 
and r = 0.4 meV, parameters that arc now established to 
be typical for MgB2. Based only on the low-temperature 
data, a conclusive decision is difficult, as both fits, though 
different in detail, capture the major features of the data 
(with the exception of the zero-bias peak [[7l| ). 

Figure ^ finally shows the temperature dependence for 
such a junction. Even at high temperature, a clearly de- 
veloped gap structure is visible in the raw data, which 
together with the obvious closing of the peak position 
along with the closing of the gap is conclusive evidence 
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FIG. 4: Temperature dependence of an SIS spectrum showing 
exclusively the small gap pl| . 

of the SIS nature of this junction. To further analyze 
such data, a standard BCS SIS model was employed, 
where F was adjusted to best fit the lowest tempera- 
ture data and then held constant with increasing tem- 
perature, leaving Ag the only parameter used to adjust 
the fits to higher temperature data. These fits closely 
reproduce the evolution of (i) the peak height, (ii) the 
zero bias conductance peak, and (iii) the filling of the 
gap caused by thermally activated quasiparticles. Thus 
all important features connected to the coherence peaks 
are captured by this model, and we may safely use it to 
determine As(T) from this data with high precision. 

The results of such fits, performed on three different 
junctions, are shown in Fig. ^. The small error bars rep- 
resent the uncertainty in fitting the data, and there is 
excellent reproducibility among the three sets of data. 
This evolution of A 5 closely resembles the theoretical 
prediction (see Fig. H) for a two-band model, and is in- 
compatible with this gap being an isolated order param- 
eter, as -^^fr would be significantly lower than the weak- 
coupling limit. This temperature dependence thus gives 
direct evidence for the presence of a second, larger order 
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FIG. 3: Low ternperature conductance spectrum (symbols) 
9] along with fits to an SIN (dashed line) 
(solid line). 
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FIG. 5: Temperature evolution As(T) from BCS fits to three 
independent sets of data |4l| ]. 
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FIG. 6: Conductance spectrum showing direct current con- 
tribution from the 2D band along with a fit to an SIS model 
using a weighted sum of two BCS DOSs. Inset shows high bias 
on an expanded y-scale. The sharp zero-bias peak reflects a 
Josephson current in this low-resistance junction. 

parameter, although typically no contribution from the 
second, 2D band is observed in SIS data. 



SPECTROSCOPY OF Al 

Such current contributions from the 2D band are ex- 
pected to show up at eV = in SIN junctions and 
at As + Al and 2Al in SIS junctions. However, it is 
necessary to very carefully analyze whatever extra fea- 
tures are observed. Both the spectra in Fig. || and Fig. 

, respectively, show additional features at bias close 
to A5 + Al, yet the shape of these features is entirely 
different. Whereas in the first case a distinct peak is ob- 
served, the second spectra displays a dip at about the 
same position. 

The extra peak in Fig. |^ can easily be attributed to di- 
rect tunneling contributions from the 2D band. The rela- 
tive weight of such contributions depends on the density 
of charge carriers in the respective bands (the 3D band 
accounts for 58% of the total DOS ||^), and the orien- 
tation of the junction. Since the small gap is hosted in 
a 3D band, it is expected to contribute to junctions of 
any given orientation, whereas the large gap requires the 
junction to be aligned with the 2D bands which generally 
will not be the case. We believe, that the latter condition 
is crucial for the observed dominance of the small gap. 

Figure ^ gives an example for the rare case that such 
alignment is achieved in SIS break junctions. Two clearly 
developed peaks are visible in the spectrum, that can be 
ascribed to 3D-3D and 3D-2D tunneling. Note again the 
qualitative difference between this second peak and the 
As + A^-dip feature that will be discussed below. 

To fit this data, we use BCS DOSs for both bands, 
adjusting the gap magnitudes to fit the well developed 
coherence peaks and choosing the mixing ratios in both 



electrodes to reproduce the observed peak heights. Such 
a model yields As = 2.3 meV and Al = 7.6 meV, con- 
sistent with our previous data (in the case of A5) as well 
as with other experimental reports. 

The fit shown in Fig. ^ uses weights of 2 and 4%, re- 
spectively, for the large gap contribution in the two elec- 
trodes . While the gap magnitudes (which contain the 
important, intrinsic information) are well-defined from 
the data, this is not the case for these mixing ratios, 
as there is no reason to assume them to be equivalent 
in both electrodes. The presence of the second peak 
gives evidence for contributions to the total current from 
the 2D band in at least one electrode, however, it does 
not give evidence as to how this 2D contribution is dis- 
tributed over both electrodes. In principle, only one elec- 
trode needs to show alignment with the 2D band to pro- 
duce a second peak at As + Al- The decisive information 
would be the presence of a third peak at 2Al, however, 
its amplitude is expected to be negligible unless the the 
2D band contributions are significant in both electrodes 
fzsl . We checked our data for such a third peak and find a 
very weak shoulder at the corresponding position that is 
consistent with 0.1% contribution to the total spectrum 
at most (a blow-up of the high-bias data is given as an 
inset to Fig. This is consistent with the weights of 2 
and 4% chosen for the fit (that yield 0.08% contribution). 
Anyway, the mixing ratio is determined by the relative 
orientation of the grains and does not contain any intrin- 
sic information. We therefore accept the aforementioned 
consistency as sufficient. 

TWO-BAND FITTING OF SIS SPECTRA 

Most tunneling work that determined A^, or at least 
found evidence for its presence, did so by studying direct 
tunneling contributions from the 2D band in SIN junc- 
tions (double-peak structures). In the previous section 
we illustrated how such contributions do, even though 
rarely, infiuence SIS junctions. However, the more fre- 
quently seen extra feature in SIS junctions is entirely 
different in shape and origin. Figure [t] gives an exam- 
ple for a well-developed dip near As -I- A^, in contrast 
to the peak at this position that we discussed before. 
This feature is commonly seen in our SIS junctions, and 
can certainly not be explained by any parallel tunneling 
channel, as such extra channels should exclusively yield 
additional conductance. 

Comparison of SIS fits and the actual data (see e.g. 
Fig. ^ reveals that (with increasing bias) the data first 
exceeds the BCS expectation (shoulder) and then distinc- 
tively drops below it (dip) fz^ , a behavior reminiscent of 
strong-coupling effects, where spectral weight is redis- 
tributed in a similar fashion. However, at this energy 
there are very few phonons available and the electron- 
phonon coupling to the 3D band at any rate is expected 
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FIG. 7: (a) Low temperature spectrum showing a clearly de- 
veloped dip structure near As -I- Al. (b) Evolution of the 
background with increasing field (0 — 1.4 T with a 0.2 T in- 
crement). The encircled As -|- A_l -structure vanishes with 
increasing field. 



to be weak. To establish this feature to be intimately 
related to superconductivity, we studied the whole back- 
ground structure in more detail. 

The bottom panel of Fig. |^ shows a magnification of 
the conductance background in the junction shown on 
top, displaying a variety of small structures, where the 
A5 -I- A^-feature by far is the most pronounced one. 
With increasing magnetic field, the superconducting gap 
structure weakens ||7^ and along with it, the Ag -I- A^- 
feature vanishes. At the same time, the other background 
structure remains unchanged, proving the high stability 
of the junction and that the A5 -|- A^-dip is closely con- 
nected to superconductivity. 

We believe, that this extra feature is a signature of 
quasiparticle transfer between the 2D and 3D bands. 
To explore this, we need to abandon the simple BCS 
model used before and employ a slightly more elaborate 
DOS, taking into account the interband quasiparticle in- 
teraction. The coupling between both bands in principle 
may be mediated by both pair and quasiparticle trans- 
fer, however, only the quasiparticle contribution yields 
corrections to the spectral shape of both energy gaps. 
Since the original treatment of two-band superconduc- 
tivity Q only considers pair transfer, it does not suffice 
for our purpose. Quasiparticle transfer was included into 
a later theoretical two-band model [p5| , resulting in extra 



0.0 



2.0 



1.5- 



^3 



1.0 



0.5: 



0.0 
0.5 



- 












: (a) 





- 














1 



1.5 



energy 



2.5 



FIG. 8: (a) Partial and total DOS in a full two-band model 
(reprinted from Ref. josj. Fig. 3, copyright 1977, with permis- 
sion from Elsevier Science), (b) Equivalent DOS calculated 
using McMillan's tunneling model of the proximity effect. 



structure in each DOS at the position of the respective 
other gap (see Fig. ||). To keep modeling mathemati- 
cally simple, we use an observation pointed out in Ref. 
|66|, viz. that the DOS in a pure BCS two-band model 
1 64 1 and the DOS in McMillan's tunneling model of the 
proximity effect are different only in that the latter 
includes both pair and quasiparticle contribution on an 
equal footing. This means, that the McMillan model may 
be used to generate the DOS in a two-band model. The 
model requires the solution of two simultaneous equa- 
tions for the respective gap functions 



l + T,/^Ai{E)^{E- 



ph 



r2A,{E)/^Ai{E)-{E~zri 



(4) 



(5) 



were are the intrinsic pairing amplitudes in both 

bands, ri_2 are scattering rates related inversely to the 
times spent in each band prior to scattering to the other, 
and r* 2 are smearing parameters in both bands which 
were added to account for lifetime effects. These gap 
functions are then used to create the DOS from the stan- 
dard BCS expression eq. (^. SIS spectra are created 
using the usual convolution eq. (|^). 
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FIG. 9: Low temperature spectrum showing direct tunneling 
exclusively from the 3D band, and an additional dip feature 
near As + Al. Thin dashed line is a fit to a (one-band) BCS 
model, solid line to a two-band model (McMillan's equations) 
The parameters used are Af" = 0, A^'' = 7.2 meV, Fi = 
4.0 meV, F2 = 1.0 meV, F^ = 0.5 meV, and F2 = 1.0 meV. 
See text for details. Inset: Derivative of the small-gapped 
DOS shown in Fig. |b. 



To further verify the vahdity of this approach, we re- 
produced the DOS calculated from a full two-band model 
|l65| using McMillan's model ||6^. Fig. ^ compares both 
results and an almost perfect agreement is found over 
the entire energy range. We therefore may safely use 
this model to fit our results and investigate the effects of 
interband quasiparticle transfer on the DOS. 

Figure ^ gives an example for a typical low- 
temperature SIS conductance spectrum, along with a fit 
to the standard-BCS model used before (dashed line) 
and a fit using McMillan's model (solid line). Whereas 
the BCS fit fails to reproduce the higher bias features, 
McMillan's model results in a close fit over the entire 
voltage range, including the dip near A5 -I- fz^ . To 
qualitatively understand the origin of this dip-feature, 
consider the explicit form of the derivative of eq. at 
zero temperature: 



dl_ 

dV 



UN 



N{E) 



eV 



dN{E + eV) 
deV 



dE+cjNN{^eV)N{{)). 

(6) 

The second term yields only a small correction (and van- 
ishes if there are no states at zero energy), thus the major 
contribution consists of a convolution of the DOS and its 



shows the derivative of the 
% displaying a pronounced 



derivative. The inset to Fig. |9 
small-gapped DOS from Fig. 
minimum just to the left of A^ (reflecting the shoulder 
in the DOS at the same position). Convoluting this with 
the DOS, that shows a sharp onset at As, yields a dip in 
the conductance slightly shifted from Ag -I- A^ towards 
lower energy, as we indeed find in our data (see also Fig. 

0a)- 

The success of this model confirms that the A5 -f A^- 
dip reflects the quasiparticle coupling in between both 
bands. To our knowledge, this so far is the only piece 



of experimental evidence from tunneling spectroscopy to 
not only show two gaps, but to also give direct evidence 
for their coupling from a detailed analysis of their spec- 
tral shape. 



RULING OUT PROXIMITY EFFECTS 

However, the applicability of this model a priori does 
not prove two-band superconductivity to be the origin of 
the dip-feature, as McMillan's model originally was de- 
veloped to describe the proximity effect, which therefore 
needs to be considered as a possible cause. In the prox- 
imity effect, a thin surface layer of reduced (or no) intrin- 
sic superconductivity is influenced by an adjacent super- 
conductor to show enhanced superconductivity (or some 
supercondutivity at all). McMillan's equations specifi- 
cally describe the proximity effect in the presence of an 
insulating barrier of thickness, d/, in between the bulk 
superconductor and a surface layer of thickness, dn- 

Allowing for such an unknown surface layer (which in 
the case of MgB2 may easily be envisioned to be e.g. 
MgOa;) the observation of a small gap is conceivable with- 
out considering two-band superconductivity. Note, that 
even the temperature dependence of A5, or more specif- 
ically its closing at Tc, is not inconsistent with the be- 
havior of a proximity sandwich. In spatially extended 
junctions (like point contacts or planar junctions) it is 
furthermore plausible to find more than one gap, as the 
surface layer may vary over the junction area (or simply 
not be present in some areas) — such spatial variations in 
sample properties over the junction area have been sug- 
gested to explain point contact data. 

There are two strong arguments against this proxim- 
ity scenario and in favor of two-band superconductiv- 
ity. Firstly, the unusual reproducibility of the small gap's 
magnitude in a wide variety of experiments using differ- 
ent samples and tunneling techniques. Referring again to 
Tab. ||, it is obvious, that almost all experiments would 
agree to As = 2.0 — 2.8 meV. To fully appreciate this 
agreement, one has to keep in mind how sensitively the 
induced gap would depend on the properties of an as- 
sumed proximity sandwich, d^ enters the induced gap 
with a mere power law, however, di, — representing a 
tunneling barrier — enters exponentially (for an experi- 
mental verification of this aspect of McMillan's model, 
see Ref. [|9j). This means, that very slight changes in 
the barrier thickness result in a variation of the induced 
gap that may easily cover several orders of magnitude. 

Again, our SIS geometry allows us to strengthen this 
point. Since the junction becomes rather complex in this 
model, we refer to the sketch in Fig. In general, the 
SIS junction is now represented by two proximity sand- 
wiches separated by the main barrier, I. Each sandwich 
in turn consists of the bulk superconductor, S, a barrier, 
I1.2, and a surface layer, Ni_2- As the surface layers form 
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FIG. 10: Schematic representation of an asymmetric SIS' 
junction caused by randomly, and non-uniformly formed sur- 
face layers. Here, the surface layers are considered part of the 
grains. See text for details. 



randomly, there is no reason to assume the 1^ and to 
be equivalent. Likewise, the effective gap values entering 
the "main" SIS junction may be different, resulting in an 
asymmetric SIS' junction. 

Such junctions are known to show a unique spectral 
shape, including a so-called "difference peak" at Ai — A2 
Jfzt . Vice versa, the absence of this feature (as seen in 
Fig. ^ gives direct evidence of the symmetry of the junc- 
tion. Testing for very small asymmetries requires a de- 
tailed analysis of the in-gap spectral shape, since the dif- 
ference peak (or zero-bias peak in symmetric junctions) 
is subject to thermal broadening — in stark contrast to 
the coherence peak (see Fig. ^). Therefore, the differ- 
ence peak may be veiled and show up only as a charac- 
teristic change of the in-gap spectral shape. Analyzing 
simulated SIS' spectra allowing for increasing difference 
in gap values, we find that 20% difference is about the 
maximum that might still be considered consistent with 

we estimate 



our data. Using the numbers from Ref. |69 



that this translates into a reproducibility of barrier thick- 
ness, I1.2, on the order of ±0.5 A. Considering the initial 
assumption that such proximity sandwiches need to form 
randomly, and independently, such a reproducibility cer- 
tainly proves this consideration to be highly unlikely. 

A second argument against the proximity effect is the 
observation of both gaps simultaneously at a single point 
in STM studies, which is very hard to reconcile with 
the assumption of variable surface layers. Although in 
principle any proximity induced, small gap structure is 
expected to show an additional feature at the intrinsic, 
large gap energy, this feature is expected to be much 
weaker than the main, induced gap structure. It cer- 
tainly is inconsistent with two peaks of equal magnitude, 
or even the high-bias peak being more pronounced than 
the low-bias peak (see e.g. data in Ref. Q). 

The appearance of two distinct gap values, the temper- 
ature dependence of the small gap, and the appearance of 



additional spectral features are therefore not only more 
easily understood in the framework of two-band super- 
conductivity, but — beyond that — not readily consistent 
with the proximity effect. 



CONCLUSION 

We have reproducibly observed a small gap, A5 — 
2.5 meV, in break-junction tunneling on MgB2 and 
traced it to high temperatures, where it closes near the 
bulk Tc- Only in rare cases, we also observe direct tun- 
neling contributions showing a large gap, A^ = 7.6 meV. 
These findings give evidence for, and are interpreted in 
terms of a two-band model. A commonly observed dip 
feature at A5 -I- A^ is analyzed using a specific two- 
band model. We argued, that this feature gives evidence 
for quasiparticle contributions to interband coupling, and 
furthermore showed it to be inconsistent with the result 
of proximity effects. Taken together with tunneling data 
from the literature, this gives very convincing evidence 
that MgB2 is one of the very rare materials showing two- 
band superconductivity. 
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